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MTIOML ADVISORT COMMITTEE FOR AEROMA.UTICS 
TECHNICAL NOTE 2200 

A SGTODY OF SECOND-ORDER SURERSOHIC-FLOW THEORY 
By Milton D. Van Dyke 


SUMMARY 


An attempt is made to develop a second, approximation to the solution 
of problems of supersonic flow which can be solved, by existing first- 
order theoiy. The method, of attack adopted, is an iteration process using 
the linearized, solution as the first step. 

I 

For plane flow it is found, that a particular solution of the iter- 
ation equation can be written d.own at once in terms of the first-order 
solution. The second-order problem is thereby reduced to an equivalent 
first-order problem and can be readily solved. At the surface of a 
single body, the solution reduces ,to the well-known result of Busemann. 
The plane case is considered in some detail insofar as it gives insist 
into the natiare of the iteration process. 

Again for axially symmetric flow the probl^ is reduced to a first- 
order problem by the ' discovery of a particular solution. For smooth 
bodies, the second-order solution can then be calcixlated by the method 
of Karmn and Moore. Bodies with corners are also treated by a slight 
modification of the method. The computing time required is several times 
that for a ceoreful first-order solution. The second-order solution for 
pressTires on cones represents a great Improvement over the linearized 
result. SecondfOrder theory also agrees well with several solxitions 
calculated by the nmerlcal method of characteristics. 

For full three-dimensional flow, only a partial particular solution 
has been found. As an example of a more genexe.1 problem, the solution is 
derived for an inclined cone. The possibility of treating other inclined 
bodies of revolution and three-dimensional wings is discussed briefly. 

INTRODUCTION 


As the linearized theory of supersonic flow approaches full develop- 
ment, the question eirises as to whether more exact approximations are 
practical.' If viscous effects are large, refinement of the perfect-fluid 
solution is useless. Otherwise, however, higher approximations are known 
to yield a closer approach to reality. In intermediate cases, an Improved 
solution is desirable in order to assess the relative effects of viscos- 
ity and nonlinearity. 



2 


MCA TN 2200 


*rhf» prototype of a hi^er-order solution for supersonic flow is 
Busemann*s series for the surface pressure in plane floj/ past a single 
hody. This single result is of considerable value in analyzing super- 
sonic airfoil sections. Two terms of the series prove sufficient for 
almost all requirements; the extension to third and fourth order is 
chiefly of academic interest. 

The atm of the present study is, therefore, to find a second 
approximation, analogous to Busemann*s result, for supersonic flow past 
■bodies which can he treated hy existing first-order theory. The natural 
method of attack, and apparently the only practical one, is hy means of 
an iteration process, taking the usual linearized result as the first 
step. Several writers have applied this procedure to subsonic flow. In 
supersonic flow, as usual, the solution is sinq>ler, so that more general 
problems cau be solved. 

This paper is a revised version of a thesis in aeronautics for the 
degree of doctor of philosophy written at the California Institute of 
Technology under a national Research Council predoctoral fellowship. It 
has been made available to the MCA for publication because of its 
general interest. 


lEERATIOn PROCEDURE 


Basic assumptions .- The problem to be considered is that of steady 
three-dimensional supersonic flow of a polytropic gas past one or more 
slender bodies. As indicated in the following diagram, the bodies are 




assvnned either to be pointed or to extend upstream Indefinitely as cylin- 
ders parallel to the free-stream direction. In either case, the origin 
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of coordinates can "be chosen so that all variations in hody shape are 
confined, to the half-space x > 0. Wind, axes are introd.uced., so that 
for X = 0 the flow is uniform and. parallel to the x-axis, with veloc- 
ity U and. Mach number M. (For definitions oj^ all symbols, see 
appendix. ) 

The bodies are slender, which means that at any point the component 
of U normal to the surface is small compared with U itself. The 
symbol e will be used throu^out as a measure of this smallness. Thus 
the ordinates of a body will be written as e times a ftinction of order 
unity. Used in this way, e serves to distinguish terms of various 
orders. 

It wl»ll be assumed that the full linearized solution to the problem 
is available. Then the aim of this investigation is to provide a second 
approximation to the exact nonlinear solution. The linearized, or first- 
order, solution is defined as the result of keeping only lineea: pertur- 
bation terms in the equation of motion. Similarly, the second-order 
solution is the result of retaining products of perturbation qtiantltles. 
In addition, however, certain of the triple products are in some cases 
found to be as important as one or more double products and are therefore 
also retained in the eqmtion. 

It may be noted that the second-order solution will not generally 
consist slnq)ly of terms of order e and though this is the case 

for plane flow. For example, the second-order solution for flow past a 
body of revolution contains terms as hi^ as e^iog^e. 

The flow is assumed to be irrotationeil and isentropic. This assump- 
tion is Justified in the first- and second-order solutions, since the 
resvilting error is found to be at most of the order of terms neglected 
elsewhere. 

Exact perturbation equation .- Under the previous assumptions, there 
exists a velocity potential fl. In Cartesian coordinates, the equation 
of motion is (reference 1, equation (39)) 

(c2 - V " 


2SlySlz^Z ~ SflgfljfUgX ” ® 0 


( 1 ) 


Here the local speed of sound c is related to Cqj its valxxe in the 
uniform stream, by 


^ co^ - + a/ + -J^) (2) 

where 7 is the adiabatic eatponent. The subscript notation is vised to 
indicate differentiation. 



h 


MCA TN 2200 


A perturljation potential $ is now introduced in the usual way. 

For convenience, however, $ is normalized throu^ division hy the free- 
stream velocity, so that 

fl = U(x + ®) (3) 

I 

The perturbation velocity at any point is then the gradient of $ multi- 
plied hy , U. 

Introducing the pertiarhation potential into the equation of motion 
gives, after some manipulation. 


*yy *zz " ^ 


XX 


y ~ ^ /pA + + + +® +* \ + 

2 X ^x y J\ XX ^yy ^zzj 

2®xx^ ®xx^^ ®yy*y^ *zz®z^ "*■ 

2®yz®y®z ^®zx®z(^ ®x) "*■ ^®xy(^ ®x)*y 




where 3 = \/m^ - 1. 

Solution hy iteration .- The exact perturbation equation (eqxiation (^)) 
is completely equivalent to the original nonlinear potential equation 
(equation (l)). Simplifying assunptions must therefore he introduced in 
order to solve it'. If if is assumed that squares suid products of the 
derivatives of $ can he neglected, the ri^t-hand side of equation (4) 
disappears, leaving the wave equation 


+ ®zz^^^ - P^^XX^^^ = ° (5) 


This equation is the basis of the linearized or first-order perturbation 
theory, so that its solution is designated hy $(!)'. 

More exact solution of equation (4) hy means of iteration was first 
suggested hy Erandtl (reference 2). The procediure has been applied to 
plane subsonic flow hy Gortler (reference 3)^ Hantsche and Wendt (refer- 
ences 4 and 5)^ Imai and Oyama (references 6 and J), and Kaplan (refer- 
ences 8 to 10). Schmieden and Kawalki (reference 11) applied it to sub- 
sonic flow past an ellipsoid of revolution. Most of these writers have 
considered the stream function rather than the potential, which restricts 
the method to plane or axially symmetric flows; The procedure is clecirly 
described hy Sauer (reference 1, p. l4o0 for the case of plane flow. 
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^ The linearized solution subject to proper boundary conditions, 

is taken as the first approximation. Substituting this known solution 
into the right-hand side of equation (4) gives 

*yy^^^ *zz^^^ - = Fi(x,y,z) (6) 

vhere Fi is a known function of the independent vsLriables. This is 
again a linear equation, the nonhomogeneous wave equation. A second- 
order solution , j subject to proper boundary conditions, can be 
sought by standard methods. The procedure can be repeated by siibsti- 
tuting 4(2) into the right-hand side of equation (4)- and solving again. 
Continuing this process yields a sequence of solutions 4(11) which, 

Tjnder proper conditions, presumably converges to the exact solution. 

This procedure bears a saiperficlal resemblance to the Picard process 
for hyperbolic eqmtions in two Independent variables (reference 12, 
p. 317). There is, however, an essential difference. In the Piceord 
process, the characteristic lines of the differential equation are known 
at the outset, since the functions F^ do not depend on the hi^est- 
order derivatives. Here, on the other hand, the characteristic surfaces 
(the Mach cones) are initially unknown. Because of the fundamental role 
played by the characteristics in the theory of hyperbolic equations ( see, 
for example, reference 13, ch. Il), it might be anticipated that 'the 
characteristics should be revised at each step of the iteration process. 
Each step but the first would then involve equations with nonconstant 
coefficients. The siibsonic countesrpart of such a procedure is known to 
converge under proper conditions (reference 12, pp. 2887289). 

However, the pixacedure outlined previously makes no provision for 
such revision. At each stage of the iteration process, the equation has 
the original characteristics of the undistiurbed flow. As a result, the 
equation has constant coefficients, which greatly facilitates solution. 
Fortimately, it will be found that this procedure nevertheless gives an 
inqproved solution neeo'ly everywhere in the flow field. 

Second-order iteration equation .- Henceforth, only the first two 
steps of the Iteration process will be considered in detail. It is 
therefore convenient to regard the second-order solution as consisting 
of the first-order solution plus a smaller additional term. Also, for 
simplicity, 4(1) will henceforth be replaced by Then 

4(2) = 0 + (p (7a) 


where 


0 = 

^ = 4(2) _ j(l) 


(7b) 
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Now 0 = is a solirtion of the homogeneoois wave equation (equa- 
tion (5))j so that substituting into the exact perturbation equation 
(equation (4)) shows that <p, as well as 4 >( 2 )^ is a, solution of the 
following second-order iteration eqmtion: 




3^(20x + 0/ + jZly2 + + J^zz) + 

+ JZ^yy^y^ + + 

2JZ5yzMz + 202 xP^2(i + JZ^) + 2JZ53^(i + j^)0y. 

Since 0 satisfies equation (5)^ the term ^0xx + 0yy + 0zz) 
the ri^t-hand side of equation (8) can be replaced by 1^0yy. and the 
eqviation for cp becomes 

1^0xx(2J^x + 0x^ + + 0z^) + 20 xx?5x + 

+ J^^yyJ^y^ + J^zz^z^ + 

20yz0y0z + 20zx0z(l + 0x) + 20xy(l + 0x)0y 


( 8 ) 


«IPyy + 9ZZ “ P^%X = 


(9) 


Here the right-hand side contains not only products of perturbation 
quantities but also triple products. The latter can be omitted for plane 
flow, since they contribute terms of smaller order (equal to those found 
in the next iteration). Otherwise, certain of the triple products should 
be retained, since their contribution is as great as that of one or more 
of the' double products and greater than any contribution from a third 
approximation. It will be seen later that triple products should be 
retained if they involve only derivatives normal to the free stream. 

Those which involve x-derivatives can apparently be neglected, so that 
the equation becomes 





- l)l^j0xxi^x ■*" 20xy^y 20zx^z 

20yz0y0z 0zz^z^ 


(10) 


Here the triple products which may be. important are grouped in the second 
line of the ri^t-hand side. 
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The adiabatic exponent 7 will be found to occur always in the 
combination 


^ (7 + 1 )M^ 

2 p2 


( 11 ) 


Making this substitution, the second-order iteration equation becomes 
finally 


'Pyy + 'Pzz - ^ 


2(N - + 


2^xy^ 2^zx0z 
+ ^zz^z^ 


+ 


( 12 ) 


Iteration eqtiation in other coordinates .- In cylindrical coordi- 
nates, eqiiation ( 12 ) becomes 




0 

2 (rr - I)p2jz5^jz5x + 4 + 


f ^xx ^‘xB 


(13) 


The terms whose fom is indicated in the last line are the triple products 
which will be fovtnd to be negligible. 


For conical flows it is convenient to introduce nonorthogonal conical 
coordinates (x,t, 0 ) where 



ilk) 


If the body itself is conical, the pertiurbation potentials are reduced 
to functions of two variables by introducing conical perturbation poten- 
tleils (reference l4) so that 


0 (x,t,e) = xj^(t, 0 ) 


' (15) 


with corresponding definitions for and cp. The derivatives are 

given by 
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= 0 - t0t 

= P^t 


^xx ^ ^tt 
l«rr = f ?tt 


^r0 = P^te 
^0x = Iq ~ 


= XC 


= 2^80 ^xr = "It ^tt 


( 16 ) 


witli the same relations connecting and and (p sind qp. 

The iteration eq\iation hecomes 

2(N - l)t20^^(^-- t^^.) - - t^t0)^0 + 

3^^tt^^+ 2 ^ 0te^t0e - ^ ^ 


(l-t2)(ptt+^ + ^=l^| 


1 

t2 


(^8 ” ■t^t0)(^ - ■t^t)^9^ t^tt^t(^ - "t^t) 


Ul7) 


Here the grouping of terms corresponds to that in equation (13). 


BoundEory conditions .- Physical considerations suggest that the flow 
should satisfy the following conditions: 

(1) The resultant velocity is tangent to the surface of the body 

(2) All flow perturbations vanish identically everywhere upstream 

of the plane x = 0 


The theory of hyperbolic differential eqtjations shows that these two 
requirements are svtfficient to determine the solution. The first iotposes- 
one condition along the timelike surface of the body, and the second 
imposes two conditions on a spacellke surface. This corresponds me^the- 
matically to the case of mixed boundaiy conditions (reference 12, p. 172) 
and leads to a determinate solution (see reference 13, p. 85 ). 

The tajagency condition may be written formally as 


Vft'VS = 0 


( 18 ) 
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where S(x,y,z) = 0 is the equation of the surface of the hody. In a 
more useful form it becomes ^ for the first- and second-order problems. 


— = Slope (on the STarface) (l9a) 


— — — = Slope (on the surface) (19^) 

1 + + <Px 


Here means the cross-wind component of the normal derivative of 0 

at the surface of the body. In plane flow or for planar systems 
is 0yj and in axially, symmetric flow 0c slope of the 

body is measured with respect to the free-stream direction. If the 
first-order tangency condition (equation (I9a)) is satisfied exactly, 
the second-order condition can be simplified to 


2s. = Slope (on the surface) (l9b‘) 

* • 

In linearized theory, the tangency condition (equation (l9a)) is 
frequently approximated by neglecting con^rison with unity. If 

the corresponding approximation is made in the second-order problem, the 
two tangency conditions become 

0c = Slope (on the surface) (20a) 

% 

^c 9c _ Slope (on the surface) (20b) 

This approximation will not be made except for plane flow, since other- 
wise it apparently causes unnecesse^Tr loss in accuracy. 

A planar system is defined to be a system for which the first-order 
tangency condition cein be imposed at a plane parallel to the free stream 
rather than on the siurface of the body (reference 15, p. 52). Thin flat 
wings are planar systems, while slender pointed bodies of revolution are 
not. For planar systems the second-order tangency condition can also be 
imposed at the plane, provided that the value of 0y is calculated at 
the surface of the body (0x and q>x may be calculated at the plane). 
That is, for plansa: systems the tangency conditions are 


(^y)plane ( Slope )^1 + ^)piane 


(21a) 
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^^^plane “ (Slope) ^1 + + %c)plane " (^V) surface (21b) 

1 

Corresponding results hold for quasi-cylindrical bodies, which eire bodies 
of revolution whose radius varies so slightly that the tangency condi- 
tions cem be inqposed at a circular cylinder parallel to the free stream. 

The other two boundary conditions are that 


0(q,y,z) + <p(0,y,z) = 0 

^x(Ojyjz) + qRx(o^y^z) = o 


( 22 ) 


These conditions must be satisfied by the first-order solution alone and 
must therefore be satisfied also by the additional second-order potential 
alone. Consequently, the conditions are 


0(O,y,z) = 0 
^(0,y,z) = 0 


(23a) 


«p(0,y,z) = 0 
q>jc(0,y,z) = 0 ^ 


(23b) 


Determination of pressure .- When the potentieJ. field has been deter- 
mined, the net velocity q at any point is given by 

q2 = (U + u)^ + v^ + (24) 

where 



( 2 ) 

X 



y 


( 2 ) 

( 2 ) ' 





( 25 ) 
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in "both Cartesian and cylindriceQ. coordinates. Because the flow was 
assumed to be isentropic, the pressure coefficient is given by 



1 + 2L_=Ji #/i _ 

I 

7-1 

- 1 


2 

. 


where p^ and hre the free-stream pressure and density. 

It is the practice in liiiearized theoiry also to linearize the pres- 
sure relation. Substituting equation (24) into equation (26) and 
eacpanding gives 


U V^ + W^ O u^ p u v^ + w^ 

n = .2 - - + m2 H ^ + 

^ U xj2 U2 U u2 



In linearized theory only the first term is ordinarily retained. This 
is satisfactory for plane flow or flow past planar systems, since the 
contribution of the remaining terms is truly of higher order. In fact, 
for plane flow past a single body it happens that the next two terms 
cancel identically. 


For slender bodies such as a cone, however, orders of magnitude 
are not so clearly distinguished. Busemann suggests (reference l4) that 
the second term is then srofficiently large compared with the first that 
it should be retained also. This view is supported by Li^thill (refer- 
ence l6), who shows that the resulting solution is correct up to the 
order of the quantities contributed by the second term. Again, the 
third term, which also involves squares of perturbation quantities, is 
conqjarable with the second at hi^ Mach numbers and might logically be 
retained. Having gone this feir, it may be simpler to use the exact 
relation. > 

Each of these four possibilities for the first-order flow past a 
5° cone is compared with the exact solution (reference 17) in figure 1. 
The series (equation (27)) is seen to alternate in this case. It con- 
verges so slowly, however, that linearizing the pressure relation intro- 
duces much greater errors than linearizing only the equation of motion. 
Adding each of the quadratic terms in turn causes fluctuations nearly 
as great as the error due directly to nonlinearity in the equation. 

The point of view to be adopted here is that calculating the veloc- 
ities and calculating the pressure are two essentially distinct 
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operations. ; A certain degree of approximation may be necessary in order 
to solve for the velocities, but the pressure relation need, not then be 
approximated, to the same extent simply for the sake of consistency. For 
it may happen that the resulting errors (thou^ of the same mathematical 
order) are greater than those d.ue to the original approximation. Indeed, 
this is evidently the case in the first-order solution for a cone and 
will be foTuad true to a greater extent in the second-order solution. 

Moreover, in the second-order solution so many terms of equation (27) 
must be retained that it is usually simpler to use the exact relation. 

For these reasons, the exact pressure equation (equation (26)) will be 
used throu^out except in the case of plane flow. 


FLAKE FLOW 


The second-order solution for conditions at a single surface in 
plane supersonic flow was given by Busemann ( references l8 and 19 ) • By 
using the iteration procedure, the solution will now be found throu^out 
the flow field, including the Case when several bodies interact. 

The solution for plane flow is of interest chiefly insofeur as it 
serves as a guide in more complicated problems. In peurticuleo:, it pro- 
vides insist into such details of the iteration process as the question 
of its success and the effect of sharp comers. 

Role of a particiilar solution ..- The second-order iteration equation 

can be attacked by standard methods, and in the case of plane flow a 
solution can be found directly. For plane and axia].l y symmetric flows, 
however, a particular solution of the iteration equation can be written 
down at once in terms of the first-order solution. This solves the 
problem, because the conq)lete solution consists of a particular solution 
plus a solution of the homogeneous equation, and the latter can be 
obtained by existing methods. That is, the additional second-order 
potential may be written as 

<p = t + X ■ (28) 

where 

i]f any particular solution of nonhomogeneous iteration equation 

X a correction potential which is a solution of correspon^ng 

homogeneous wave equation and serves to correct the tangency 
condition 

The problem for X is the usual first-order problem whose solution is 
assumed to be known. 
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The role of the particular solution is to transfer the nonhomogeneity 
in the problem from the equation, where it is troublesome, to the boundary 
conditions, where it can be handled by existing theory. For linear par- 
tial differential equations it is always possible in principle to trans- 
fer nonhomogeneities in this way from the equation to the boundary condi- 
tions, and vice versa, by adding a suitable function to the dependent 
variable (see reference 20, p. 236). 

Since the particular solution ilf will be found in terms of the 
first-order solution, it will vanish upstream of the pleuae x = 0. Then 
the correction potential must also vanish there, so that two boundary 
conditions are given by 


X(0,y,z) = )^(0,y,z) = 0 (29) 

From equation (l9b')j the tangency condition for X is found to be 

ijf- + Xp 

— = Slope (on the surface) (30) 


or, in the case of planar systems, from equation (21b) 

(’^y * ’^)plane ' ^ (31) 

It should be noted that, ailthou^ q> is small conq)ared with 0, 
this is not necessarily true of either i|f or X alone. 

Particxilar solution for plane flow .- The first-order equation for 
plane flow is 

^yy - P^i^^xx = 0 (32) 

The general solution is 

0(x,y) = H(x - 3y) + j(x + Py) (33) 


where H and J are functions chosen so as to satisfy the boundary 
conditions. 
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In the iteration eguationj all triple products cem he neglected, 
and eqiiatlon (12) becomes 

*1^ “ ~ ^xy^J (3^) 


It can he verified diregtly that a particvilar solution of this equation 
is given hy 

Ilf = - ly + 1 (35) 

To this must he added a solution X of the homogeneous equation (equa- 
tion ( 32 )), which hah the form 

X = h(x - py) + J(x + py) (36) 

where h . 6Uid J are functions determined hy the second-order boundary 
conditions. 

For flow past a single boundary (such as one surface of an airfoil) 
the first-order potential (equation (33)) contains only one or the other 
of the functions H eind J. In this case ^xy^y “ P^^xx^x that the 
iteration equation reduces to 

cpyy - p2<p^ = 2J^p2ujzl^j^^ (37) 


The particular solution may then he simplified to 


Ilf = m 2 I yjzj^jzly ‘ (38) 

and the correction potential contains only h or J, according as the 
first-order solution contains only H ’ or J. 

Flow past a cTirved wall . - As an example of the application of the 
particular solution, consider flow past a wall which at some point 
begins to deviate sli^tly from a plane (see the following figure). 

The wall can he represented by 


y = eg(x) 


(39) 
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where c is a parameter small compared with unity and g(x) is a 
continuous function of order unity which vanishes for x s 0. 


y 


, I 

I 




This is a planar system,, so that the tangency conditions are given 
hy equation (21), The approximation of equation (20) can also he made. 
Consequently, the first-order pi*ohlem is 


f*yy - = 0 

0y(X,O) = €g»(x) 


m 


0(o,y) = = 0 

The solution is 


0 = -|- g(x - Py) (4l) 

p 

Substituting into the ri^t-hahd side of equation (34) gives the 
iteration equation 

qpyy - = 2M2jre2g>(x - Py)g*’(x - Py) (42) 
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According to equations (38) and (36), the solution is 


cp = t + X = e2y[g*(x - py)] ^ + h(-x - py)' 




Imposing the approximate second-order planar tangency condition (equa- 
tion (2Ib)) gives 


h’(x) = |g* (^ )] ^ - g(x)g”(x) > 


(44) 


so that 


h(x) = 


g(x)g'(x) + I [g'(l)]® M 

2p2 

lI 0 


(45) 


The cQnq)lete second-order perturbation potential is therefore 


$(2) = _ py^ _ g2j 


g(x - py)g*(x - Py) + ^ y|g’(x - py)]^ 


M^(H - 2) 
2p2 



[g‘U)]2 di 


(46) 


Th6 same result can he found hy solving equation {U2) directly, using 
the tnq>ulse method (reference 12, p, l64). 


On the surface of the wall, the streamwise velocity perturhation 
is given hy 

The pressure coefficient at the wall can now he calculated from equa- 
tion (27) which, upon replacing N hy its value from equation (U), 
gives 

Cp = I eg'(x) + ^ ~ [Gg'(3c)]2 (1*8) 

This is the well-known result of Busemann (references l8 and 19). To 
second order, the surface pressure coefficient depends only upon the 
local slope. 
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Role of cliaracterlstics , - It was pointed out previously that, 

because of the -underlying significance of the characteristic surfaces 
for solutions of hyperbolic equations, it might be expected that the 
characteristics would have to be revised successively at each stage of 
the iteration. However, an iteration process -was chosen which permits 
no such revision. It is therefore pertinent to inquire in this sinq)le 
exan^le what role has been played by the original and the revised 
characteristics . 

Only one of the two families of characteristics -will be considered. 
The original characteristics of this family are the lines of slope 


dy _ 1 

dx 3 




These are the Mach lines of the undisturbed flow which run downstream 
from the -wall (see the preceding diagram). They epre also characteristics 
of equation (32) in the mathematical sense (reference 12, ch. 5; refer- 
ence 13, ch. II). 


It can readily be shown that, if the first-order streamwise perttirba- 
tion velocity at any point in a flow is. u^^J, then the revised local 
valioes of Mach ni^er and p are given to first order by 


m( 1) = M 1 + p2(N - 1) 


U _ 


(50a) 


p(i) = \/m( 1)^ - 1 = p 


3 


1 + m2(N - 



(50b) 


By using this result together -with the first-order solution (equa- 
tion (4l)), the revised downstream Mach lines are found' to have the slope 


dy _ 1 
dx p 


1 + 6g.(x - py) 

P 


(51) 


These are not the mathematical characteristics of the iteration equa- 
tion (equation (^2)) for the reason that fractions of the hi^est-order 
derivatives have there been transferred to the rl^t-hand side and 
regarded as known. Mathematically, the characteristics continue to be 
given by equation (^)., 

Physically, the characteristics are lines along which discontinu- 
ities in veloci-ty deri-vatives are propagated, and this definition is 
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completely equivalent to the mathematical one (reference 12 ^ p. 297 ) • 
Therefore, in the second-order solution derived above, discontinuities 
in acceleration must occur along the original characteristics. 

Suppose, however, that no such discontinuities occur. For flow 
past a single body the downstream characteristics are also lines along 
which the velocity is constant, provided that shock waves do not appear. 
Setting 


^ *xy^^^ dy = 0 


d®y(2) = ^ + *yy^^^ dy = 0 


it is seen that the velocity is constant if 


( 52 ) 




( 53 ) 


For the second approximation (equation ( 46 )) the velocity is constant 
along lines of slope 




( 54 ) 


which, according to equation (50b), are the revised characteristics. 
Consequently, althou^ the characteristics have not been revised in the 
mathematical sense, the solution behaves physically as if they had, so 
long as discontinuities do not occvur. The question of discontinuities 
will be considered in the next section. 


The connection between the original and revised characteristics 
can be interpreted physically. The right-hand side of the Iteration 
equation may be regarded- as representing the effect of a known distri- 
bution of supersonic sources throu^out the flow field. The influence 
of this soTirce distribution spreads .downstream along both f a milies of 
original characteristics. The resulting velocity changes are just such 
that the second-order velocities become constant along the revised 
rather than the original characteristics. 

Finally, it is Interesting to note that the second-order potential 
is constant on lines which bisect the original and revised character- 
istics. For, setting 




dy = 0 


( 55 ) 
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$(2) foiand to "be constant along lines of slope 

£ = ( 56 ) 

Flow past a corner and a parabolic bend . - A simple case in wMch 

discontinuities may ocean: is that of flow past a sharp corner. The exact 
solution is known to involve an oblique shock wave with attendant veloc- 
ity discontinuities for cos^resslon and a continuous Prandtl-lfeyer fan 
for expansion. 

Denoting the tangent of the deflection angle by e, positive for 
compression (see the following figure), the function g(x) appearing 


y 



in equation (39) is 




( 57 ) 


From equation ( 46 ) the second-order perturbation potential is fooind to 
be 


®(^^(x,y) 


-i. (x - Py) ■+ ^ (x - ^y) 

P P*^ 


w2ir 


e^x 


( 58 ) 
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to the ri^t of the. line x = 3y and zero to the left. Consequently, 
in either compression ( e > O) or expansion ( € < O) the second-order 
potential .suffers a discontinuous drop along the Mach line from the 
comer, of strength proportional to the distance from the comer. Such 
a discontinuity cannot he admitted, which indicates that the iteration . 
process fails in this region. 

In the case of ccm^ression, the solution can he corrected hy 
analytically continuing the perturbation potential upstream until it can 
he joined continuously to the free-stream potential. (This is permis- 
sible since the line of discontinuity is not actually a characteristic.) 
From the result of equation (56) the juncture is seen to occur along the 
line from the comer which bisects the tq)stream and downstream Mach 
directions, as indicated in the following diagram: 



Mach lines hefore-and after adjustment 
of potential discontin\iity. 


The adjusted discontinuity corresponds to a shock wave, for it is known 
that em oblique shock bisects the Mach lines to first order (refer- 
ence 13 , p. 35^). In the case of eaqjansion, this type of correction 
cannot be 'justified, since it would involve continuation of- the free- 
stream potential across a true characteristic. Instead, a Prandtl -Meyer 
fEui must be inserted. 

Evidently the Iteration process is successful except within an 
angular region of order £ lying near the Mach line from the comer. 

In particular, the pressvire is given correctly everywhere on the surface 
of the wall. 
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It Is enll^tenlug to olaserve that the alternative method of iter- 
ation, in vhich the characteristics are successively revised, fails in 
the same region. The potential is doiible-valued over a fan-shaped region 
in the case of compression and is left undefined over a similar region 
in the case of ea^jansion (see the following diagram). The same artificial 



corrections are necessary to coi^lete the solution. 

Consider next flow past a paraholic hend which is represented hy 

y = i ex2 X ^ 0 ■ (59) 


Prom equation ( 46 ) the second-order perturbation potential is found to 
be 


*^^^(x,y) = (x-Py)2- 5Ll5±|I:il g2(x_ py)2_ ^ g2y(x_ py)2 (6o) 


The potential and also the velocities are continuous, so that the 
previous difficulties do not occur. The acceleration is discontinuous 
across the original characteristic x = Py, which in this case happens 
to be also a revised characteristic. However, a new complication arises. 
It is well-known that, in the exact solution for the compressive cage. 
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the characteristics form an envelope, as shown in the accompanying figure. 



First- and second-order flow past a parabolic bend. 


Inside the cusp the potential is triple-valued (reference 13, p. Ill), 
so that a shock wave must be inserted. This envelope must also arise in 
the second approximation, since the characteristics are no longer 
parallel. However, the second-order potential given by equation (6o) 
is single-valued, so that it cannot predict the formation of an envelope. 
Again the iteration process fails in a part of the flow field. 

It can be seen that the alternative iteration process, using revised 
characteristics, will produce an envelope. 

Convergence for plane flow .- The examples just considered indicate 
that the success of the Iteration procedure should be carefully investi- 
gated. A step of an iteration process may be considered successful if, 
in some sense, it significantly improves the solution. In particular, 
one is interested in the success of the second-order solution. 

It should be noted that a divergent process may be successful for 
many steps and that, on the other hand, convergence does not necessarily 
imply success. In practice, however, one would expect a convergent 
process to be successful. As used here, success is a subjective ndtion, 
not amenable to analysis. Consequently, only the convergence of tte 
iteration procedure can be considered in any detail. 

Unfortunately, proofs of svtfflclent conditions for convergence have 
not been obtained, even in the case of plane flow. However, the above 
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examples suggest certain conjectures regarding convergence. These will 
be stated and some arguments for their plausibility will be advanced. 

For flow past a slightly curved plane wall represented by y = eg(x) 
the solution obtained by iteration using the revised characteristics is 
conjected to converge in any bounded region adjacent to the wall provided 
that 


(a) € is sufficiently small 

(b) g(x) is continuously differentiable 

If g(x) has only a piecewise continuous derivative, the convergence 
holds except possibly in fan-shaped regions springing from each comer, 
which lie near the original Mach line and subtend an angle of order e. 


For the Iteration process actvially adopted, in which the character- 
istics are not revised, the first n steps are conjectured to form peirt 
of a convergent process provided that 

(a) 6 is sufficiently small 


(b’) g(x) has continuous derivatives up to (n. - l)st order if 
the potential is required, nth order if the velocities are required 

If condition (b’) is satisfied only piecewise, the result holds except 
possibly in fan-shaped regions springing from each comer. 


In the first case, condition (a) is necessary in order to insure 
that the solution be vinique, as is clear from the example of the para- 
bolic wall. The above examples also show that condition (b) is necessary. 


If the sufficiency of these two conditions is assumed, their connec- 
tion with condition (b’) in the second case can be illustrated by analogy 
with a mathematical model (suggested by Dr. C. R. DePrlma) which retains 
the essential difference between the two iteration processes - namely, 
that the correct characteristics are not used in the method actually 
adopted. Consider the first-order problem given by equation (^); 

^yy " ^xx ~ ^ 1 


J^y(x,0) =€g’(x) ► 

0(o,y) = 0x(o^y) = 0 


(61) 


where P = 1 has been taken for convenience. The solution (equation (4l)) 
was 
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0 = -€g(x - y) 


(62) 


Now It is attempted to solve this prohlem using characteristics which 
differ from the true characteristics hy 0(e). Thus consider the 
equivalent problem 


0yy - (1 - O0X2C = 

« 

0y(X,O) = €g*(x) •- 


( 63 ) 


0(O,y) = 0x(O>y) = 0 


and solve hy iteration. In the first approximation the right-hand side 
can he neglected^ so that 

0yy - (1 - €)0^ = 0 (64) 

which has the solution^ subject to the boundary conditions^ 

= -6g(x - \/l - ey) (65) 

Substituting this into the right-hand side of equation (63) gives the 
iteration equation for the second approximation: 

^yy - ” ®^*^xx = -€^g”(x - Vi - cy) (66) , 

\ 

Using the impulse method (reference 12, p. l64) gives the solution, 
siibject to the boundary conditions, 

= -eg(x - yjl - €y) + ^^yg'(x_- \/l - ey) (6?) 

But this is Just the Taylor series expansion, correct to 0(e2)^ of the 
true solution (equation (62)). Subsequent iterations add additional terms 
to the expansion. Hence, despite the use of slightly incorrect charac- 
teristics, the iteration process converges to the correct solution. 

The connection between conditions (b) and (b’) is thus seen to be that 
the existence of sufficiently many continuous derivatives compensates 
for the fact that the wrong characteristics are used. 
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AXIALLY SYMMETRIC FLOW 


Before discussing the general solution for "bodies of revolution, it 
is convenient to consider the simple problem of a cone. In this case 
the second.-order solution can he found, directly. The results will he 
useful in indicating which triple products should he retained in the 
general case. 

Flow past a cone .- Consider flow past a slender cone of semivertex 
angle tan'lgj as shown in the following diagram: 



Flow past a cone. 


The flow is conical and eocially symmetric, so that the iteration equa- 
tion is given hy equation (17) with ©-derivatives omitted. Including 
the boundary conditions from equations (I9a) and (23a), the first-order 
problem is 

(l.t2)0^t+i=O 

= €[l + - Be^(B€)][ 


m 
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Using "the Integrating factor t/ - t^, the equation can he Integrated 
to give the well -known result that 

^ = -A^sech"^ - \/l - t^) (69) 

which Is understood to vanish except within the downstream Mach cone, 
where t = 1. The tangency condition Is satisfied by putting 


e2 

\/l - p2g2 + 62sech"l(3€) 


(70) 


At the Mach cone (t = l), all velocity perturbations vanish, so that 
the first-order solution predicts no deflection of the shock wave from 
the Mach cone (see reference 13, p. ^3). 

Substituting the first-order solution Into the Iteration equation 
(equation ( 17 )) gives 


(l - = a 2# |2 i + 2(H - 1) •- p2A (Tla) 

and from eqmtlons (l9h’) and (23b) the corresponding boundary conditions 
are 

p^(p€) = «[^(pe) - p6^(pe)] (71b) 


qj(oo) = (p^.(oo) = 0 


(71c) 


Equation (71a) can again be solved using the Integrating factor t//l - t^. 
The various Integrals encoxintered can Inveurlably be treated* by Integrating 
by parts one or more times. . Using the second boundary condition, the 
complete conical second-order perturbation potentlsuL Is fovuad to be 


■^2)(t) = -A^sech"^ - \jl - t^)+ A%4^ B^sech"^ - \Jl - t^) + 


3 / 2 " 

(sech"^)^ - (H + 1) Vl - ^2 sech-^ - ^ ~ 



(72) 
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From equation (l6), the strearawise and radial velocity perturbations are 
^ = -A sech”^ + A^l^ |b sech"^ + (sech“^)^ - (N - l) - 

U L \/2~^ 

(N + 1) - ^ ^£^^1 (73a) 


g = A i£^ t A%2 Pb - g i£I F secylt , ^ + 1) 1 + 

P U *b L. w ^ 

(H . 1) t_se^ ^ 1 (a ^ 

VirnE2 t3 j 


The constant B must be adjusted so as to satisfy the tangency condition, 
equation (71b). In actual computation it is easier to adjust B numeri- 
cally in exactly this fashion rather than to calculate it from the cumber- 
some e:q)ression which could be written down. The pressure coefficient at 
any point can then be calculated from equation ( 26 ). 

The laat.term in the bracket in equation (71a) is the triple 
product which is retained in iteration equation (17)» Its 

retention ie now Justified by noting that its contribution - the last 
tenn in equation (72) - is of the same order as other terms near the 
’surface of the cone (t = pc). Actually, it contributed a second term, 
which has been neglected since it is at most of order c^sech'^pe. It 
can also be verified that the other triple products, whose form is indi- 
cated at the end of equation ( 17 ) ^ are in fact negligible, since they 
contribute at most terms of order €^( sech~Ipc)2. Consideration of a 
further iteration ^dicates that a third approximation would add terms 
no greater than e°(eech"Ipe)3, which is greater than the terms Just 
neglected. 

The second-order resvilt for surface pressure coefficient is com- 
pared in figure 2 with the exact solution (reference 17) for cones of 
5°, 10°, 15°^ and 20° semivertex angles. Also shown for comparison are 
the first-order resiilts based upon the exact e^qiression (equation ( 26 )) 
for the pressure coefficient. The second-order solution is seen to pro- 
vide a much better approximation throughout the range of Mach numbers up 
to the point at which the Mach angle equals the cone angle, beyond which 
the perturbation solutions have no physlceLL meaning. 
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Shock-vave angle .- The solution for plane flow past a comer 
suggests that the seconcL-order solution for the cone may fail neeir the 
Mach cone. However, if it is valid there, a first approximation to the 
shock-wave deflection and, consequently,, the entropy change can he calcu- 
lated from the fact that to first order an oblique shock bisects the Mach 
•lines. It ■was noted before that first-order theory predicts no differ- 
ence bet'ween the shock position and the }&ch cone. 


Assume provisionally that the solution is valid at the Mach cone, 
while indicating by ? the possibility that it is not. From equa- 
tion (73) the velocity perturbations jiist behind the Mach cone are 


/'uN I _2M%6^ 
vuyt=i 

('I') I 23M%6^ 

'‘U/b=l 


(74) 


so that the perturbation is normal to the Mach cone. Here A (equa- 
tion ( 70 )) has been -approximated by €^. Prom equation (50b) the cotan- 
gent of the revised Mach angle Just behind the cone is found to be 

3|l - 2M^(U - 1)6^] (75) 


The upward stream inclination there is approximately 
the Mach lines have the slope 


^ , so that 

\uA=i' 


g I i(l . ■ (T6) 


Therefore, the slope of the shock wave differs from that of the original 
Mach cone by 


1 ? M^^ ^ (7 + 1)^M^ 

" P P 435 


(77) 


This problem has been treated rigorously by LigdJthill (reference 21) and 
by Broderick (reference 22), who find that actually 


tan X 


1=3 (y + 1)^ g4 

P 8 p5 


(78) 


which is 


ll times the above 
2 


result. 


The discrepancy- means that the 
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second-order solution fails near tbe Mach cone. It seems renarkahle 
that the result is in error only to the extent of a consteurt factor. 

The entropy increase through a weak oblique shock wave is propor- 
tional to the cube of its inclination away from the Ifech lines. Conse- 
quently, the entropy rise through the shock wave from a cone is O(e^), 
as noted by Lighthill (reference 21). 

Particular solution for axially synnnet'ric flow .- Consider flow past 
a body of revolution which is either a slender pointed body with nose at 
the origin or one which extends indefinitely upstream with constant 
radius a for x = 0 (see diagram). Q?he latter shape corresponds to 


r 



Flow past bodies of revolution. 


the external flow past a sharp-edged, open-nosed body with supersonic 
internal flow. With slight modification the subsequent development can 
be applied to internal flow as well. The meridian curve can be repre- 
sented in the first case by 


r = R(x) = eg(x) 


X ^ 0 


(79a) 
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and in the second hy 


r = R(x) 



cg(x) 


X ^ 0 
X ^ 0 


(79t) 


Here e is again a parameter aman compared with unity, and g(x) is 
a continuons function of order unity which vanishes for x = 0. 

The first-order prohlem is 



( 80 a) 

^(x,R) = R‘(x)[l + 5^x(2c,R)] 

( 80 b) 

^(0,r) = ^x(0,r) = 0 

( 80 c) 


The solution is known to he (reference 23) 


5ZJ(x,r) 



F(l) dg 



g)2 > p2r2 


ncosh"^ 


x-h 

Pr 

F(x - Pr cosh u) du 


Jo 


( 81 ) 


The second form is \iseful for carrying out differentiation, after which 
the first form can he restored. The derivatives which will he required 
are 


0x 



x-h 

PF- 

F'(x - Pr cosh u) 


du 



F»(l) 


( 82 a) 



X-h 

pr 

F'(x- pr cosh u) cosh u 


du = 



:-Pr 


(x- l)F'(l) dl 

V(x-|)2.p2r2 


( 82 b) 
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With coordinates as shown in the preceding diagram, the lower limit of 
integration h is 0 for the iwinted hody and -3a for the semi- 
infinite "body. The function F(x) may he interpreted physically as the 
strength of a supersonic line source along the x-sixis. It is determined 
hy the tangency condition, eqiiation ( 80 b), which gives the following 
integral eqiiation of Volterra type for F’ : 


px-PR(x) 

nx-pR(x) 

(x- 5 )P'(S)d! 

^ I F»(l) dg 

V(x- 1)2- p2b2(x) 

[ \Z(x-g)2-p2B2(3,) 

Jb ^ 

Jb 


(83) 


be 


The second-order iteration equation is 


found from equation (l3) to 


<Prr + ^ - P^'Pxx = 


2(N - + 20xrlzJr + 

‘J^rr4^ + 

®(^xx^x^^ J^xx^r^^ ^xr^x^r) 


(84) 


The solution for the cone suggests that the terms indicated in the last 
line are negligible. 

It will now be shown that a particular solution of this equation 
is given by 

Tjr(x,r) = + Nr^) - ^ l#r^3 (85) 


The first group of terms contributes the first line in equation (84), as 
can be verified by direct substitution. The last term in equation ( 85 ) 
accounts for the term 0rr0r^ follows: 
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r I; ■ " "J^(| + i + i r^rr^r + 

^ ^ — I P2r0xr^ - I 32r0xxr0r) 

= JZ5r 0rr?^r-|rjZ5r(0rrr + %^-^-32iZ5^ - 
+ I r0rr)(0rr + ^ " P^i^xx) “ 

P^(i ^xxJ^rr + ^xJr ' | ^i^xr^)] 

= 0rr0r^ ~ r^xx^ + ^xx^r - 1 ( 86) 


vhere repeated use is made of the fact that 0 satisfies equation (80a). 

The last group of terms consists of triple products involving x-derivatlves, 
which have already heen neglected in equation (84), so that the result is 
proved. 

The correction potential X' is a solution of equation (80) and can 
he written as 


nx-pr 


X(x,r) = - 


f(l) d| 


X.-1 3C-h 

pcosh”-L ^ 


Jb 


V(x - 5)2 - p2r2 


f(x - Pr cosh u) du (87) 


vJO 


Using equation (82) the second-order tangency condition (equation (l9h*)) 
is found to he 


iIf;p(x,R) 


/ 


nx-PR 

— 

nx-PR 

(x - e)f’( 1) d5 

i|rx(x,R) - 

' f*(0 d| 

V(^- 5)^- 

Jb 

}J(x- 1)2- p 2 r 2 
"b 


( 88 ) 


which is again a Volterra integral eqtation. 

Methods of solving integreJ. eqaiation . - Discovery of a particular 
integral for bodies of revolution reduces the second-order problem to 
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the same form as the first-order problem - namely, the solution of a 
Volterra integral equation. Various methods of attacking this problem 
are listed by Hayes (reference 15, p. 1^0). 

■ An indirect method consists in assuming that the -unknovn source 
strengths in equations (8l) and (87) can be represented by a few terms 
of a polynomial, for exanrple, that in eqmtion (81) 

F(x) = + . . . + (89) 

The resulting solutions were introduced in a more formal manner by Hayes 
(reference 15, p. 38) ^ who has discussed their properties in detail. The 
first term alone gives the -potential for the cone, equation (69). Addi- 
tional terms give the solution for simple families of shapes. However, 
the method is not suitable for bodies having discontinuities in slope or 
curvature. Consequently, a more direct procedure is desirable. 

Karman first introduced an asynqrfcotic solu-t^Lon of the integral equa- 
tion (equation (83)) which has come to be known as the slender-body 
approximation. For slender bodies, the source strength F(x) appearing 
in eqviation (8I) is found to be approximately proportional to the rate 
of change of cross-sectional area. Thus 

F(x) » ^ = R(x)R*(x) (90) 

where S(x) = nR2(x) is the cross-sectional area of the body. Li^thill 
has shown (reference 16) that if R(x) and its first two derivatives are 
of order e, eind R’ is continuous, then this determination of F(x) is 
correct to the order of terms retained in the first-order solution. For 
purposes of the second-order solution, it can be shown that F(x) may be 
determined in this way only if R” is also continuous. This means, that 
the body must have continuous curvature, which is a severe limitation. 
Moreover, the slender-body approximation is foimd generally to cause 
unnecessary loss of accxiracy even though the mathematical order estimate 
of the error is small. Consequently, this approximation is not to be 
recommended if it can be avoided. 

The most satisfactory way of solving the integral equations is to 
use a step-by-step numerical procedure. In first-order theory the usual 
method, introduced by Keamian and Ifoore (reference 24), is to assume that 
the unknown source distribution can be approximated by a polygonal graph. 
This is equivalent to super in^o sing a number of conical source lines of 
different strengths, each shifted downstream with respect to its prede- 
cessor, as Indicated in the following diagram: 
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Equivalence of polygonal source and svmi of conical sources. 


The latter viewpoint is more convenient for computation. The strengths 
of the source lines eire determined in succession hy satisfying the 
tangency condition at a series of points on the surface of the body. 

The details of this procedure are cleeirly explained in reference 1, 
page 77. 

For piorposes of a second-order solution, this procedure must he 
modified in one respect. The source distrihution F(x) must not he 
approximated hy a polygon unless it dctually has comers. The reason 
for this is that a comer corresponds locally to adding a conical source 
line, which gives the solution for a cone. But it was fomd in the case 
of the cone that to second order the velocities are discontinuous across 
the Mach wave. This would cause false pressure Jump in the flow field. 

Instead, the procedure must he carried out iislng source lines of 
quadratic strength. The source strength F(x) can then he approximated 
smoothly, so that false pressure Jvmrps do not occur. A single source 
line of this type represents the flow past a slender pointed hody with 
a cusped nose (see sketch), as is clear from the slender-hody approxi- 
mation, equation (90). 



> Body formed hy source line of quadratic strength. 
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Method of solution for smooth. Tjodles .- The second-order solution 
will he described first for bodies having continuous slope. Modifica- 
tions for treating shair comers will be discussed in the next section. 

The procedure is indicated in the following diagram. The axis is 




Ifethod of solution for smooth bodies. 

divided into inteivals by choosing points with abscissas at each of 

which a source line is to begin. Good accuracy is usiially obtained if 
the inte^al length is not greater than 3 times the local radius. The 
tangency condition will be imposed on the surface of the body at the 
points Pn^ which lie on the Mach lines from the points at l^. 

For pointed bodies, the first-order solution is started with a 
conical source from the origin which gives the proper conical tip. This 
potential and the derivatives which are reqiilred are 
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0o = -Cox(sech“H - \/l - 

(^o)x ~ “Co secli~^ 

(0o)r = PCo 

u 

vhere 


(i^o)xx 

(0o)xr 

(0o)rr 


_ Cq 1 
X ^ - t2 

_ 3Cq 1 
P^Co 

* t2C^t^ 


Co = 


\/l - p2g2 + 62secll“l^e 


(91) 


and the semivertex angle of the conical tip is tan"le. No such term is 
required for the semi-infinite body. 

The subsequent procedure is the same for either body. Quadratic 
source lines are started from each of the points Ij, and so forth. 

Por the pointed body is also at the origin, while for the semi- 
infinite body it is at -pa. For the nth such source line, the poten- 
tial and its derivatives are given by 

0n = rCn(x - 5n)2|j^l + i. Tn^^sech'^n - 

(0n)x = -2Cn(x - Sn)(sech”^n - ^1 - 

(?fn)r = PCn(x - sech'^n^ 


(0n)xx = “2 Cb sech“^n 


(0n)xr = 2 PCb 




(^)rr = “P^Cnl 


- Tj 
\ Tn2 


+ sech”^i 


(92) 



where 
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The cons'tan'ts are determined successively hy imposing the first- 

order tangency condition in turn at each of the points Pn. From equa- 
tion (80b), the condition is that 

n-1 

^ (^^)r 

— = Slope (93) 

1 + X Wx 

where the sinnmation begins with n = 0 for the pointed body and n = i 
for the semi-infinite body. In this way, values of the complete first- 
order potential 0 and its first and second derivatives are calculated 
at each of the points Pq. 

The velocities due to the particular second-order solution’ ilr- cein 
then be calculated at the same points. Differentiating equation (85) 
gives 

fe = - I 

• 

Finally, the second-order correction potential X is determined by 
repeating the procedure used for finding new constants such that the 
second-order tangency condition is satisfied. From equation (l9b’), the 
condition is that 

n-1 

■’l^r (^n)r 

Slope (95) 

(^)x 

The second derivatives of X need not be calculated. 

The complete second-order perturbation velocities are the sums of 
the. contributions from 0, i|r, and X. Then the pressure coefficient 
be calculated at each point P^ from equation (26). 

The con^Juting time required is several times that for a careful 
first-order solution. 



can 
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Treatment of “bocLles with corners .- Suppose that the meridian curve 
of the hody has a sharp comer, which, for convenience may he assmned to 
lie on the Mach cone from the origin ( see fig. ) . Then the method of 


r 



solution must he modified for two reasons. 

In the first place, the intervals between source lines would have 
to he chosen extremely smun in order to obtain an accurate first-order 
solution behind the comer. This difficulty can he eliminated, however, , 
by adding a new solution which causes a sharp deflection of the stream- 
lines. In this way the comer is effectively removed. 


Such a solution can he found by approximating to equation ( 83 ) in 
the vicinity of the Mach cone. The resulting Abel integral equation can 
be solved to show that, in general, a potential having discontinuous 
nth derivatives resvilts from a source distribution along the axis which 


is initially proportional to 
tion (81) gives 


n-J 

ir ^ . 


Setting F(x) = x 


n-| 


in equa- 


px-Pr 


0(x,r) = - 


n-| 

i ^ dl 


(x - pr^^ 


r, 1 

(1 - l)“-5 at 


\/(i - 1)2 - a2r2 




1 + 

2Pr 


) 


( 96 ) 
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This integral represents the analytical continuation of the hyper- 
geometric function, so that, except for a constant factor. 


jZl(x,r) = (x - F^|, |, n + 1, - 


X - Br' 


where a is the radius at the comer. The potential is understood to 
vanish except within the downstream Mach cone from the origin. The 
hypergeometric .functions occtirring here can all he expressed in terms of 
complete elliptic integrals with real modulus. 

QThe solution for a comer is obtained by taking n = 1. Then 

, P T 1(98) 

; I # - =(rj4)] 

^rr ‘4i§'T^ "M)] 


Here t is the conical variable introduced in equation (l^), and 

■'(Hr) are the complete elliptic integrals of the 

first and second kinds wlth.modxilus k such that k^ = i— 

1 + t 

From the tangency condition, equation (80b), it can be shown that 
the above solution should be multiplied by 

(R2* - Rl')[l + (^x)J /gc 

P + R2’ . " 


4o 
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in order "to cEuicel "the corner. Here Ri’ and' R 2 ’ are tlie slopes of 
the meridian curve just "before and after the corner and ($x)l 
value before the comer. The first-order solution can then be continued 
as described in the previous section. , 

The second difficulty which arises is that, . even though the first- 
order solution may be exact, the second-order solution described previ- 
ously is incorrect behind the comer. It is clear that the local pres- 
sure jump at a comer should have the value correspon din g to plane flow 
past a comer of the same angle. However, the second-order solution 
described above yields the first-order, rather than the secOnd-order,- 
value of the pressure jimip. The method of solution must be modified in 
order to obtain the correct result. 

The proper procedure would be to fibd the solution for the case when 
the corner has been slightly romded off and then pass to the limit of a 
sharp comer. However, the following simpler procedure is found to give 
exactly the same result. 

The particular solution ijr calculated from equation ( 85 ) is discon- 
tinuous along the Mach wave springing from the comer. If the disconti- 
nuity vanished at the comer, the solution could subsequently be revised 
as In the case of plane flow (see diagrams of Mach lines in section 
"Flow past a corner and a parabolic bend"). However, there is a finite 
jimp in ijr directly at the comer, which cannot be allowed. Conse- 
quently, the correction potential X must involve an eqxxal and opposite 
jump. A potential having stich a discontinxiity is obtained by setting 
n = 0 in equation (97)* Then 



Adding a suitable multiple of this potential cancels the discontinuity 
in ilf. The second-order solution can then be continued as described in 
the preceding section. It can be verified that the pressure jump at 
thd comer has then the correct second-order value. 

It is instmctive to analyze the behavior of a comer from another 
viewpoint. It was pointed out before that the ri^t-hand side of the 
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Iteration equation may be considered to represent the effects of a known 
distribution of sources throti^out the flow field. In the case of a 
sli^tly rounded corner, this source distribution will be weak except 
between the Mach lines from the comer. As the corner shrinks to a 
point, the source intensity will increase in tl^t region in such a way 
that the total strength remains constant. In the limit, the source 
distribution will behave like a Dirac 6-function along the Mach line 
from the corner. The particular solution for plane flow (equation (35)) 
takes account of this Inqmlsive function, so that the coirect solution 
is automatically obtained. In the case of axially symmetric flow, 
however, it is clear that the particular solution given by equation (85) 
misses the contribution of the impulse. It is therefore necessary to 
correct this shortcoming by adding the step-function potential given by 
equation (lOO). 

Coi!g)arison with numerical solutions .- The accTiracy of the second- 

order solution for bodies of revolution can be evaluated by con5>arison 
with examples ceilculated using the numerical method of characteristics. 

The first body to be considered is a 12.5-callber ogive, which has 
a semivertex angle of 16.26° at the tip. The second-order solution was . 
calculated for this body at a Mach number of 3.2^. This represents a 
severe test of the method, because the Mach angle is then only 10 percent 
greater than the tip cone angle. Intervals were chosen such that the 
points Pn lay at 0.1, 0,25# 0.5# 1# 2# and 3.5 calibers (fig, 3)» The 
pressure distributions calculated by first- and second-order theory are 
compared in figure 3 with the results of various characteristics solu- 
tions sianmarized in reference 25* The second-order solution appeurently 
coincides with the characteristics solutions to within the accxiracy of 
the characteristics method. 

The second body to be considered consists of a cone of 10° semivertex 
angle followed by a circvilar cylinder. The characteristics solution foy 
this body at a Mach number of 2,075 has been given by Liepmann and Lapin 
in reference 26. The first- and second-order solutions were calculated 
beyond the comer using the modifications discussed in the preceding 
section. Figure 4 shows the shape of the body, the location of source 
lines, and the velocity distributions calculated by first-order theory, 
second -order theory, and the method of cliaracteristics. Again, the 
second-order results agree with the characteristics solution to within 
the accuracy of the latter method. 

Series expansion with respect to thickness .- An alternative method 
of solving the exact perturbation equation (equation (4)) by successive 
approximations is to assume that the solution can be expanded in powers 
of the thickness peirameter e. Thus the exact perttirbation potential 
is written as 
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k2 


(1) (2) _ (3) o L 

$= $£ + 9 e + ®eJ+ 

Substituting into equation (U) and equating llie powers of 6 
seqTience of equations 


( 102 ) 
yield a 


( 1 ) ( 1 ) Jl) 

*yy *zz - p *xx - 0 


( 2 ) ( 2 ) ( 2 ) 

®yy *zz “ P ®xx 


of 

2t^ (N - l)p2 


( 1 ) ( 1 ) 

*xy®y "■ 


( 1 ) ( 1 ) 


zx 




U (103) 


which can he solved in succession. The first is again the usual linear- 
ized equation. This method was applied to plane subsonic flow in refer- 
ences ij- and 9. 

Schmieden and Kawalki first pointed out (reference 11 ) that the 
power series assumed here does not always exist. In general, terms of 
the form e™(log €)“ begin .to appear in the third-order solution for 
plane flow and in the second-order solution for three-dimensional flow. 
This difficulty can be met by assuming a more general series of the form 


® = 


( 1 ) 

® e + 


( 2 ) , 
9 


(3) . 

$ e- 


® e^og 


(5), 

e + ® e 


(104) 


On the basis of this assunqrtion, Broderick has developed a' secoiid- 
order solution for supersonic flow past slender pointed bodies of revolu- 
tion (reference 27). The eumlysis is rather lengthy, since the simpli- 
fication resulting from the discoveiy of a particular solution does not 
appear. The results are limited to shapes for which the cross-sectional 
area is given .by an analytic function, or at any rate possesses contlnuoiis 
derivatives up to the fourth order. This is a severe limitation since, 
for example, the two bodies discussed in the previous section are not 
admissible. 

Broderick’s result can.be obtained by expanding the present second- 
order solution in powers of t and log t for small valvies of. t. The 
logarithmic terms arise from the series 
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sech.-lt = log I - ^ . (0 < t < 1) (105) 


The expansion will now he carried out for the case of flow past a cone. 

It is clear from equation (70) that the constant A in the first- 
order solution (equation (69)) is given approximately hy 


A = + . . . 


(106) 


Substituting this value into equations (73a) and (73b), expanding in 
powers of t and log eind imposing the tangency condition, equa- 
tion (Tib), shows that 


B = - 


2m2 - 




(107) 


Then according to equation (73), the velocity perturbations on the sur- 
face of the cone are 

. . (108a) 


Z = € _ e^iog — + 

U e xog + . . . 


(108b) 


Replacing R by its value from equation (ll), eqijation (27) gives for 
the pressure coefficient on the surface of the cone 


a(a log ^ - 1) + 

:''|3p2(log . (5^^ 

isy] 




(109) 


This is Broderick's result (reference 27, equation (81)). 
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This series is compeirecL with the original form of the second-order 
solution in figure 5- For the most slender cone, the expansion in series 
causes only a moderate loss in accuracy. For more practical thicknesses, 
however, the expansion reduces the accuracy to such an extent that for 
the cone of 20° semivertex angle, Broderick’s polutlon is inferior to the 
first-order result. The reason must he that the -iteration process itself 
converges more rapidly than do the subsequent expansions which are 
required to reduce it to series form. Hence, terminating all expansions 
at terms of the order of those retained in the iteration process results 
in an unnecessary loss of acctrracy. 


THREE-DIMEasrSIOmL FLOW 


Part->ftT particular solution .- It mlgpit he hoped that a particular 
sol\rtion, which so greatly sinpllfies the iteration for plane and axially 
symmetric flows, could he found for the general three-dimensional case. 
The -various methods of existing first-order theory could then he applied ' 
Immediately to the problems of second-order flow past such shapes as 
inclined bodies of revolution and three-dimensional -wings. 

A part of such a particular solution is found at once, being common 
to the -fcwo special cases. Consider the three-dimensional itera-tion 
equation (equation (12)): 


«Pyy + «Pzz “ P^%CX = ^ 


2(,H - l)^^|Zlxx^Z5x + + 20zx0z 

?^yy^^ + 20yz0y0z + 0 zz0z^ 


(no) 


It can he readily verified that taking N = 0 and neglecting the triple 
products in the last line, a particular solution is given by 

r = ( 111 ) 


which appears in both equations (35) and ( 85 ). 

The iteration equation is thereby reduced to 


^ + ^zz- P^<5Pxx ^ 


23%|ZW!x + 

0yy^^ + 20yz0y^Z + 0 ZZ^Z^ 


( 112 ) 
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It has not "been possible to find a particular solution of this equation 
in terms of the first-order potential. The solutions for plane and 
axially symmetric flow do not appeso: to suggest a generalization. On 
the other hand^ there is no assurance that such a solution cannot he 
found, so that one is tempted to search fvirther. If the triple products 
are neglected, the ri^t-hand side of equation (112) vanishes for 
7 = -1 (N = O). However, investigation of the previous solutions indi- 
cates that the idea of here taking 7 = -1 is not legitimate. 

• 

In the absence of a complete peurbiciQar integreJ., the reduced iter- 
ation equation (equation (112)-) must "he attacked by more conventional 
methods. In principle, it is always possible to find a particular solu- 
tion of a linear nonhomogeneous equation with the aid of the fundamental 
solution associated with thd differential operator. For the three- 
dimensional wave operator which occotb here, the fundamental, solution is 


■■ H :.. ■ (113) 

\/(x - 1)2 - p2|jy _ 1^)2 + (z _ 


which can be interpreted as the potential at any point (x,y,z) lying 
Inside the downstream Mach' cone from a unit supersonic source at 
With the aid of Green's foirmula, it can be shown that a particular solu- 
tion of 


^ + <Pzz P^xx “ ®*U,y,z) 


illh) 


is given by 


^Kx,y,z) 


2^ 

2ic 


FU,ti,^) d| dq dt 

V(x - t)2 - p2[(y - n)2 + (z - 5)2] 


(115) 

ft 


where the integration extends throu^out that ’portion of the forward 
Mach cone from the point (x,y,z) within which F is defined. 

In practice, the integration indicated in equation (115) is gener-^ 
ally not feasible. For exanq>le, even the simplification of axial symme- 
try reduces equation (U5) only to a double integral of F(x,r) multi- 
plied by a complete elliptic integral of conqpllcated argument. Avoiding 
such integrals by disdovery of the particular solution clearly represents 
a great simplification in this case,^ 

IComparing the two methods would lead to the evaluation of definite 
integrals involving complete elliptic integrals, which mi^t be of some 
Interest. 
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In the following sections, one example of a three-dimensional solu- 
tion will he given, and the possibility of treating other shapes will he 
discussed thereafter. 

TnrUnftd cone .- The problem of a cone at an angle of attack illus- 
trates the use of separation of variables to reduce the three-dimensional 
iteration egiaation to tractable form. 

Two alternative coordinate systems are suitable for bodies of revo- 
lution at an angle of attack. In wind axes the body is inclined, while 
in body axes the stream in^jinges on the body obliquely. The latter 
system is simpler for first-order problems and is probably better for 
the second-order solution also. However, wind soces will be used here, 
since otherwise the iteration equations must be rederived. 

To facilitate in^josing the'tangency condition, it is convenient to 
apply an oblique transformation (see, for example, reference 28 , p. l8). 
This effectively unyaws the axis of the body (but distorts the surface) 
while leaving the wave operator vinchAnged. Thus three • different 



Coordinate systems for. inclined cones, 
coordinate systems are required, as shown in the above sketch: 
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Wind axes: 

Body axes: 

Oblique axes: X,'Y,Z X,E ,0 X,T ,0 

the latter two being used also in cylindrical and conical form. 

To simplify the solution, it will be assumed that the angle of 
attack a is so small that its square can be neglected. This will give 
a solution nonlinear in the body thickness but linear in a and will 
therefore yield the correct initial slope of the lift curve. Then the 
three systems of coordinates are related according to the following table: 


z = ^ = Z 


I = X - ay 


Tj = y + ax 




X = 5 + ctlAi 


Y = Ti 


To this approximation 


R = p 


(116a) 


= .(i- 


Ttt COS dj 


(116b) 


the surface of the cone is 


(ll6c) 


Ts = 


(U 7 a) 


Tg = - 3#e2a cos •3 


( 117 b) 


and the velocity components are related by 


® - T$^ 


(118a) 


$p = - T®j}a cos 0 


(118b) 
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^ ^ - T$Qi)a sin © (H8c) 

where, as in equation (15)^ the conical potential is introduced hy 

(119a) 
(1191)) 

first-order problem, referred to oblique coordinates, is found 

to be 

(l-T2)fc. + ^ + |§=0 (120a) 

3^ + cos 0 + a cos 0 = + ^ - T^) at T = Ts (l20b) 

0(«,0) = ^(*,0) = 0 (120c) 


4(X,R,e) =X*(T,0) 

where 



The solution is the sum of potentials for a conical line source and 
dipole (reference 1, p, 7^) and has the form 



Vl - t2^ + _ T 


sech"^ a 


cos 0 


(121) 


Substituting into the tangency condition (equation (I20b)) and expressing 
values of fxinctions on the cone in terms of their values at T = pe by 
TtiPHTiR of Q?aylor e3q>ansions, it is found that 


A = 


\jl - p^^ + e^sech"^! 


1 + J^A 


C = pc2 




p2^ 


pe 


- sech 




(l + 2e^)\Jl - p2e2 + p2e2gech“^pe 


(I22a)‘ 


(122b) 



The streamwise perturhatibn velocity is 


02 = (0 - T^t) + P 0 ra cos 0 


+ (2C + BA) ^ 


= -A sech“^ + (2C + PA) — — a cos 0 

T 


(123) 


Then according to equation (ill), the partial particular solution is, in 
conical form. 


= I m 200 z = m 2 a|a sech"^(sech"^ - \jl - T^) - 

t 

f(sech"^)^Ji 


(3C + PA) ~ _ (2C + PA) ~ 


CT( sech"-^)*" |a cos 0 
= ^(I)+ cos 0 




There remains to solve the reduced iteration equation given by 
equation ( 112 ) which, after transformation of coordinates, becomes 
« . 

(1 - t 2)^ + ^ ^ = i^|asr^^(0 - T^) + 

PT^( 3 T^ - 20)a cos 0J + p20j^,^ \ ( 125 ) 

Substituting equation (l2l) into the ri^t-hand side gives 


T 


(1 - t2)^ + ^ ^ = Am2Ja[ 2N - p2A )- 


Vl - t2 


2U(C + pA) i? - + IT(2C + pa) i - 2p2AC 

T\/l - t2 ^ t5 


a cos 0 / ( 126 ) 


2 
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This is reduced to two total differential equations by setting 

■ ^(T,0) = ^^^(T) + cos 0 (127) 


Therefore the complete second-order potential consists of a term indepenr 
dent of a plus one proportional to a cos 0. The first of these must 
he the previous solution for the symmetrical cone (equation- (72)), so 
that ^I) is known. The equation for is 

(l - t 2)^(II) + = -2Al#r2N(C + pA) + 

^ ^ I T\/l - t2 


H(2C + pA) ^ - 2p2AC 


(128) 


Setting 

^II)(T) = T£d(T) (129) 

reduces this to a linear first-order equation in to which can he inte- 
grated to find that 


^II)(T) = 


- T sech 


(C + PA)HT(sech-^)^ + |> P^AC 


+ (3C + 


2PA)H ^ + 


(130) 


The tangency condition -(equation (l9h’)) 
conditions 



separates • into the two 



at T 


pa 

(131a) 
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at T = p6 (I31b) 


The first of these is the previous relation (equation (Tllj)) which deter- 
mined the constant B in equation (72). Similarly, the second of these 
determines the constant D in equation (130). 

Series eacpansions for pressure and normal force .- Numerical results 

have been calculated only for the case in which the solution is expanded 
in powers of T and log Carrying out the expansion, the’constemt D 
is found to he 

Then calculating the velocity components from eqiM,tion (ll8) and the 
pressiare from equation (27) gives, on the siirface of the cone, 

a COB d + . . . (133) 

Here (Cp)o value for zero angle of attack, given hy equa- 

tion (IC^). Integrating gives the normal-force coefficient, based on 
cross-sectional area: 

This result has been obtained also by Lighthill (reference 29), who has 
calculated the lift on bodies of revolution having analytic meridian 
curves by assuming a series expemsion for the velocity potential. 

Stone (reference 30) has developed a solution for inclined cones 
which is linearized with respect to a, but otherwise exact. Kbpal 
(reference 31) has published tables of the numerical results of Stone’s 
theory. A comparison of equation (13^) with this exact theory and with 


Cp = (Cp)o - 1 - ^ | jj2 + 
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Tsien^s first-order solution (reference 32) is shown in figure 6 for 5° 
and 10° cones. The earlier discussion of series expansions suggesi^s 
that the agreement mi^t improve if the solution were not expanded in 
series. 

Shock-wave position .- If the solution were valid at the Mach cone, 
the velocity components there would he, from equations ( 72 ), (. 118 ), ( iP k) 
and ( 130 ) ; 

9 

-(§) = f (§) = ^ - 3(3N - l)3a cos d] ( 135 ) 

For sinpliclty, using equation (122), A and C have here heen approxi- 
mated hy 6 2 and pe2. Conparlng equations (7^) and (77) it is seen 
that the difference between the shock-wave angle and the Mach emgle 
would he 

* ’ 


X - sin"l i = _ 1 p%%( 3 N - l)e\t cos d ( 136 ) 


Hence the ratio of the emgular rotation of the shock wave to that of the 
cone would he 


I 1 I p2i^H(3N - 1)6^ (137a) 


It was seen previously that althou^ the solution does not in fact ' 
converge at the Mach cone, the shock-wave deflection calculated in this 
way is correct for the unpitched cone except for a factor of 1^. It 

ml^t he supposed that the same correction factor would apply here. 

Kopal (reference 3l) tabulates . val\ies of 6 /a calculated from Stone's 
theory, and from these it appears that a factor of 3 raliher than 1 ^ is 

required, so that actiially 

1 = 1 ^%%(3K - 1)6^ . (137b) 

Figure 7 shows a conpairlson of this modified result with the exact values 
for a 5 ° cone. 

It must he enqphaslzed that eqxiation (l37b) represents nothing more 
than a conjecture. It could probably he verified, however, hy extending 
the solution of Li^thill (reference 21) or Broderick (reference 22) to 
the case of angle of attack. 
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Future Investigation .- Two large classes of problems which have only 

been touched upon deserve further study. One of these is wingsj the 
other, bodies of revolution at an angle of attack. The example of the 
inclined cone was undoubtedly made awkward by the use of wind coordinates. 
The iteration equation should be rederived in body coordinates and the 
solution extended to general bodies of revolution. It is possible that 
in this form a particTilar integral could be discovered. That there is 
good possibility of success with this problem is STjggested by the fact 
that Li^thill was able to obtain a general solution by assuming a 
series expansion (reference 29). 

The possibility of discovering particular integrals of the iteration 
equation might be investigated more systematically. If none can be found 
for general three-dimensional flow, special cases such as conical flow 
should be studied. 

Possible treatment of wings .- Possibly the most useful application 
of first-order theory is to thin flat wings. Ho attempt has so far been 
made to find the second-order solution for a wing. It seems likely, 
however, that solutions can be found at least for conical problems. In 
this case the iteration equation can be reduced, by the standard conical 
theory (references l4 and 28), to the problem of solving Poisson* s equa- 
tion inside a circle. 

Two difficulties can be emticipated. First, if the wing has sub- 
sonic edges, infinite velocities arise there, so that the assunrption of 
small perturbations is violated. It is known that in first-order theory 
this is no essential objection, since the pressure is found correctly 
except in the Immediate neighborhood of the singularity, and the inte- 
grated val\ies of lift and moment are correct to first order. Kaplan 
(reference 10 ) and Schmieden and Kawalki (reference 11) have indicated 
that this result extends to the second approximation for subsonic* flow, 
so that probably no real difflctilty exists. 

Secondly, if the wing has supersonic edges, the failure of the 
iteration process along Mach lines from the apex can be expected to 
affect the surface press\ires. Again it is possible that integrated 
values will be correct to second order. Otherwise, it may be possible 
to adjust the solution in those regions, in a manner similar to that 
shown in the diagram of Mach lines in the section "Flow past a corner 
and a parabolic bend." ^ 

Higher approximations .- It seems unlikely that a third or higher 
approximation would ever be Jtistified. Other neglected factors, chiefly 
viscosity and heat conduction, should certainly be considered first. 
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However^ the Busemaim second-order result has been ejctended to third and 
even fourth order (reference 33) ^ and various writers have considered 
the third approxi mat ion for plane subsonic *flow (references 1 > «^id 8). 
If a third approxhnatlon should be considered worth while, the Iteration 
could be repeated. Again the cases of flow past a cxirved wall and a 
cone would serve as helpful examples. 

Application to subsonic flow .- The iteration eqiiatlon and the 
particular Integrals are in no way restricted to supersonic flow. The 
partic\ilar solution for plane flow might profitably be compared with the 
subsonic solutions of references 4 to 10. 

The particular solution for axially symmetric flow makes possible 
a second-order solution for bodies of revolution at subsonic speed. In 
this case, the integral equation can be treated by the methods used for 
the airship problem. 


California Institute of Technology 

Pasadena, Calif., December 19^9 
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A,B,C,D 

c 

Cn 

Cp 

E(k2) 

f(x),F(x) 

Pn(x,y,z) 

g(x) 

J J 

K(k2) 

M 

2^'^ 

P 

Pn 

q 
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APPBHDIX 

SYMBOLS 


constant reference radius for "body of revolution 

alascissa at vhich source distribution for body of 
revolution begins 

constants determined by boxmdary conditions 
local speed of sound 
constant coefficients of series 
pressure coefficient 

complete elliptic integral of the second kind with 
modulus k 

source-distribution functions for body of revolution 

known right-hand side of nth-order iteration equation 

continuous function of order unity which vanishes 
for X 0 

arbitreiry functions of one variable 

complete elliptic integral of the first kind with 
modulais k 

free -stream Mach number 


local static pressure 

points on body of revolution at which tangency condition 
is Imposed 

local speed of flow 

radius in cylindrical coordinates 
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R(x) radius of meridian curve of "body of revolution 

S(x) .cross-sectional area of tody of revolution 

t conical variatle 



u,v,w 


perturtatlon velocity con^ionents in Cearteslan or cylindrical 
coordinates 


U 


x,y,z 


X,Y,Z' 

XjR;6 

X,T,0_ 

a 

p = \/5^ 

7 

6 

e 

d 

\ 

i ' 

p 

T 

«P 

0 


free-stream velocity 

Cartesian coordinates with, x in free-stream direction 

ohllgue axes (see diagram Of coordinate systems in section 
entitled "Inclined cone") 

angle of attack 


adiahatlc exponent 

• angular rotation of shock wave on cone due to angle of 
attack 

parameter small compared with unity 
azimuthal varlahle in cylindrical coordinates 
angle of shock wave from free-stream direction 
Cartesian coordinates of variable point 

body coordinates (see diagram of coordinate systems in 
section entitled "Inclined cone") 

local density 

conical variable referred to x = | rather than x = 0 

additional second-order perturbation potential 

first-order, (lineeurized) pertiurbation potential, 
as $(l) 


same 
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$ 


$ 


(n) 


(n) 

$ 


X 


i 

lift 

ft 

03 

( 

( 

( )(“) 
n . 

Subscripts; 

c 


o 

s 

1 , 2 


exact perturbation potential 
ntb-order pertTorbation potential 

nth tem in series expansion of perturbation potential 

second-order correction potentieQ. ' 

particular solution of second-order iteration equation 

partial particular solution for three -dimensional flow 

complete velocity potential 

auxiliary variable (see equation (129)) 

result of nth iteration 

Independent of 0 (see equation (12?)) 

proportional to a cos 0 (see equation (12?)) 

conical potential; for exanq)le, * = x$ 

differentiation in cross-stream direction - component of 
normal direction which is perpendicular to free stream 

free-stream conditions 

surface of cone 

values ahead of and behind a corner 



58 


MCA TH 2200 


REPEREINCES 


1. Satier^ Rolaert: Introduction to Theoretical Gas Dynamics. Edwards 

Brothers, Inc. (Ann Arhor), 19^7. 

2 . Prandtl , L. ; AUgemeine Uberlegungen Tiber die StrcSmung zusammen- 

druokharer Fliissig^eiten. Atti dei Convegni 5, R. Accad. d* Italia, 
1936 , pp. 169-197. (Also printed in Z.f.a.M.M., Bd. 16 , Heft 3, 
1936 , pp. 129 -l^j available in translation as MCA QM 805, 1936, 
and as R.T.P. Tremslation Ho. 1872 , British Ministry of Aircraft 
Production. ) 

3 . Gortler, H. : Gasstromungen mit Dbergang von Unterschall- zu 

Uberschallgeschwindi^eiten. Z.f.a.M.M. , Bd. 20, Heft 5, Oct. 19^, 
pp. 25 ^ 262 . 

4. Hantzsche, W. , and Wendt, H. : Der KompressibilitatseinfluSs fur 

dunnfe wenig gekrOramte Profile bei TTntersc ha 1 1 Rnbt rt nfl 1 gV ei t. , 
Z.f.a.M.M. , Bd. 22, Heft 2, April 1942, pp. 72-86. 

5 . Hantzsche, W. : Die Prandtl-Glauertsche Haherung als Grundlage fur 

ein Iterationsverfahren zur Berechnung kompressibler Unterschall- 
stromungen. Z.f.a.M.M., Bd. 23t Heft 4, Aug. 1943, pp. I 85 -I 99 . 

6 . Dnai, Isao: Two-Dimensional Aerofoil Theory for Con 5 )ressible Fluids. 

Rep. No. 294 (vol. XXI, no. 9)^ Aero. Res. Inst., Tokyo ImperieQ. 
Univ. , May 1944, pp. 283 - 331 . (in Japeuiese with English abstract.) 

7 . Imai, Isao, and Oyama, Seiichi: Kie Third Approximation of the Thin- 

Wing-Eapansion Method for Co^ressible Fluids. Rep., Inst. Sci. 
and Technol. , Tokyo Univ., vol. 2, nos. 3-4, March-April 19^, 
pp. 33-44. (In Japeinese. ) 

8 . Kaplan, Carl: The Flow of a Compressible Fluid past a Curved Surface. 

MCA Rep. 768, 1943 . 

9 . Kaplan, Carl: The Flow of a ■Compressible Fluid past a Circular Arc 

Profile. MCA Rep. 794, 1944. 

10. Kapla n , Qarl: Effect of Compressibility at Hi^ Subsonic Velocities 

on the Lifting Force Acting on an Elliptic Cylinder. MCA Rep. 834 , 
1946. 

11. Schmieden, C. , and Kawalki, K. H. : Beltrage zum TtostromTmgsproblem 

bei hohen Geschwindlgkeiten. Bericht SI 3 , Teil 1, L.G.L. , 1942, 
pp. 40-68. (Available in translation as MCA QM 1233, 1949. ) 



MCA TN 2200 


59 


12. Courant, R. , and Hilbert, D. : Methoden der mathematiscben Physik. 

Vol. II. Julius Springer (Berlin), 1937. (Reprinted, Interscience 
Publishers, Inc. (Hew York), ,19^3. ) • 

13. Courant, R. , and Friedrichs, K. 0.: Supersonic Flow and Shock Waves. 

Interscience Publishers, Inc. (New York), 19^. 

14. Busemann, A. : Inf initesimale kegelige ihjerschallstroinung. Deutsche 

Akademie der Luftfarhtforschung, Dec. 4, 19^2, pp. 455-^70. (Also 
printed in Jahrb. 19^3 deutschen Akademie Luftfahrtforschung, 

Bd. pp. 105-122; available in translation as MCA TM 1100, 

1947. ) 

15. Hayes, Wallace D. : Linearized Supersonic Flow. Ph. D. Thesis, 

C.I.T., 1947. 

16. Lighthill, M. J. : Supersonic Flow past Bodies of Revolution. 

R. &M. No. 2003, British A. R.C., 1945. 

17. Staff of the Con^niting Section (tmder the direction of Zdenek Kopal) : 

Tables of Supersonic Flow around Cones. Tech. Rep. No. 1, Center 
of Analysis, M.I.T., 1947. 

18. Busemann, A. , and Walchner, 0. : Prof ileigenschaften bei TJberschall- 

geschwindigkeit. Forsch. Geb, Ing.-Wes., Bd.' 4, Heft 2, March- 
April 1933j pp. 87-92. (Available in translation as R.T.P. 
Translation No. l'^6, British Ministry of Aircraft Production. ) 

19. Busemann, A. : Aerodynamischer Auftrleb bei Uberschallgeschwindigkeit. 

Atti del Convegnl 5^ R. Accad. d* Italia, 1936, pp. 328-360. (Also 
printed in Luftfahrtforschung, Bd. 12, Nr. 6, Oct. 3^ 1935^ 
pp. 210-220; available in translation as British Air Ministry 
Translation No. 2844.) 

20. Courant, R. , and Hilbert, D. ; Methoden der mathematischen Physik. 

Vol. I. Julius Springer (Berlin), 1931. (Reprinted, Interscience 
Publishers, Inc. (Hew York), 1943.) 

21. Lighthill, M. J. : The Position of the Shock-Wave in Certain Aero- 

dynamic Problems. • Quart. Jour. Mech. and Appl. Math., vol. I, 
pt. 3 , Sept. 1948, pp. 309-318. 

22. Broderick, J. B. ; Supersonic Plow past a Semi -Infinite Cone. Quart. 

Jour. Mech. and Appl, Math,, vol. II, pt. 1, March 1949^ 

pp. 121-128. 

23. Von Karm^, Th. : The Problem of Resisteuice in Compressible Fluids. 

Atti del Convegni 5j R. Accad. d' Italia, 1936, pp. 222-277. 


6d 


MCA TBf 2200 


2k. Von Karman, Theodor, and Moore, Norton B. ; Resistance of Slender 
Bodies Mbving vith Supersonic Velocities, with Special Reference 
to Projectiles. Trans. A.S.M.E. , vol. 5^, no. 23j Dec. 15^ J-932, 

PP. 303-310. 

25. Isenherg, J. S. : The Method of Characteristics in Compressible Plow. 

Part IB - Numerical Examples. TR No. F-TR-1173C-ND, Air Materiel 
Command, U. S. Air Force, Dec. 19^7. 

26. Liepmann, Hans Wolfgang, and Lapin, Ellis; Simntiary of Character- 

istics Method for Steady State Supersonic Plows. Rep. No. SM-133^3j 
Douglas Aircraft Co., Inc., March 3j 19^. 

27. Broderick, J. B. : Supersonic Flow round Pointed Bodies of Revolu- 

tion. Quart. Jour. Mech. and Appl. Math., vol. II, pt. 1, 

March 19^, PP. 98-120. 

2&. .Lagerstrom, P. A. : Linearized Supersonic Theory of Conical Wings. 

MCA TN 1685, 1950. 

29. T.lght hm , M. J. : Supersonic Flow past -Slender Pointed Bodies of 

Revolution at Yaw. Quart. Jotir. Mech. and Appl. Math., vol. I, 
pt. 1, March 19^j PP. 76-89. 

30. Stone, A. H. : On Supersonic Flow past a Sligditly Yawing Cone. Jour, 

Math, and Phys., vol. XXVII, no. 1, April 19^^ PP. 67-8I. 

31. Staff of the Con^iuting Section (under the direction of Zdenek Kopal): 

Tables of Supersonic Flow arovuad Yawing Cones. Tech. Rep. No. 3> 
Center of Analysis, M.I.T. , 19^7. 

f 

32. Tsien, Hs\ie-Shen; Supersonic Flow over an Inclined Body of Revo- 
lution. Jour. Aero. Sci., vol. 5^ uo. 12, Oct. 1938, pp. 480-l<83.’ 

. Donov, A. : A Plane Wing with Sharp Edges in a Supersonic Stream. 

Izvestia Akad. Nauk U.S.S.R. Serie Math., 1939. 


33 















(b) 10° cone. 
Flgmre 2 .” Continued. 
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(c) 15° cone. 
Figure 2.- Continued. 
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(d) 20° cone. 
Figure 2.- .Concluded. 
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Figure 3.- Pressure distribution on a smooth body of revolution. 
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(b) 10° cone. 




0 ,- 


ConcLnued. 


MCA TN 2200 



70 


MCA TN 2200 



M 


(c) 15° cone. 
Figure 5.- Continued, 
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(d) 20° cone. 
Figure 5.- Concluded. 



(a) 5° semlvertex angle. 


(b) 10° aemivertex angle. 


Figure 6.- Comparison of varloiis solutions for normal force on a cons. 
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Figure 7.- Ratio of shock-wave rotation to angle of attack for a 5° cone. 
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